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Calculations for infinite nuclear matter with realistic nucleon-nucleon interactions sug- 
gest that the isoscalar effective mass (lEM) of a nucleon at the saturation density equals 
m* /m ~ 0.8 it 0.1, at variance with empirical data on the nuclear level density in fi- 
nite nuclei which are consistent with m* /m fa 1. This contradicting results might be 
reconciled by enriching the radial dependence of lEM. In this work four new terms are 
introduced into the Skyrme-force inspired local energy-density functional: r(Vp)^, i"^, 
and rAp. The aim is to investigate how they influence the radial dependence of lEM 
and, in turn, the single-particle spectra. 

1. Introduction 

The nuclear energy-density functional (EDF) is considered nowadays as one of the 
most promising theoretical tools to describe static properties of the atomic nuclei. 
Tremendous effort is undertaken recently to develop high-precision spectroscopic 
quality functionals. A standard way to construct the nuclear EDF is to start with 
either the finite-range Gogny ^ or the zero-range Skyrme ^ effective interaction 
and average it with the density matrix within the Hartree-Fock (HE) method. Such 
functionals may be further enriched by adding new terms. For instance, Carlsson et 
al. ^ considered a systematic generalization of the Skyrme functional by introducing 
terms up to sixth order in derivatives of the density matrix. 

Calculations for infinite nuclear matter suggest that the isoscalar effiectiwmass 
(lEM) at the saturation density should be of order ofm*/m ~ 0.8±0.lSEEIZIHl. Qn 
the other hand, description of the single-particle (s.p.) level density in finite nuclei 
requires the lEM close to unity. This contradicting conditions might be fulfilled 
together by assuming the lEM smaller than one inside the nucleus and peaked at 
the surface. Such a concept was first explored by Ma and Wambach^^in a non-self- 
consistent model and by Farine et aL^l^ within a self-consistent framework. Since it 
is impossible to obtain a surface-peaked lEM within the standard form of Skyrme 
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functional, in the present paper we examine four simple additional terms that can 
modify the radial dependence of lEM. 

The new terms are introduced in Sec. [2] Sec. [3] reports on the resulting radial 
dependence of the lEM. In Sec. [4l influence of the new terms on the s.p. spectra is 
presented. The paper is concluded in Sect. [5] In this exploratory work, we limit our 
calculations to three spherically symmetric isoscalar nuclei, ^"Ca, ^^Ni, and ^"°Sn. 



2. Extensions to the Skyrme energy-density functional 

The Skyrme energy density, H(r), consists of the kinetic and interaction parts, 

^ t=o,i 

where 

Ht = [po]p? + C^'ptApt + Clptn + C/J? + c7'ptV ■ J*. (2) 

Index t = 0, 1 denotes the isospin, and Ct are the coupling constants, of which 
one, Cf[po] — Cf + Cfpa"", depends on the isoscalar density. The potential energy 
terms are bilinear forms of the time-even densities, p, r, J, and their derivatives. The 
density 3t denotes the vector part of the spin-current tensor, 3t^\ = ^\^^v^t,^lly■ 
Readers interested in details are referred, e.g., to Ref. ^1^. 

The central field, W(r), is defined as a variation of the energy density with 
respect to the particle density, 

Z^*(r) = ffi^ = 2CnPo]pt + 2C^Ap,-fC;r,-HC7-^V.J, + <5,o ^ ^p?,. (3) 

The effective mass is proportional to the inverse of the mass field, which is a vari- 
ation of the energy density with respect to the kinetic density, 

. . 5[H{r)] h\ , r<r r ^ u\ 
——- = Mt{r) = - —dto + Clptir). 4 

The lEM is defined as rriQ/m. It is easily seen from Eq. ^ that in the standard 
Skyrme functional, the lEM is a monotonic function of the isoscalar particle density, 
Po, and therefore cannot be peaked at the surface. The term proportional to po in 
Eq. ([4]) can only decrease or increase the lEM in the interior of the nucleus. 

Thus, in order to enrich the radial dependence of the lEM, one has to add to the 
functional new terms depending on the isoscalar kinetic density, tq. We consider 
here four such terms, 

w(■^)(r) = Co^(^'')^o(r)(Vpo(r))^ (5) 
nf\r)^Cl^^r,ir)'j^, (6) 

nf\v)^Cl\T,{v))\ (7) 
H(^)(r)=C„^^''ro(r)Apo(r). (8) 
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These terms will be treated independently and dubbed as variants A, B, C, and D 
of our model, respectively. Their contributions to the mass field read, respectively, 

M[^\v) = Cl^''^'\vp,{v))\ (9) 

A^rW-C;-^, (10) 
Mf\v)^2CfT,{v), (11) 
M''f\v) = Cl'^P^p,{v). (12) 

All of these terms but Ti!"^^ contribute to the central field as well, 

W^^)(r) - -2Cl^'"'^' (Vro • Vpo + roApo) , (13) 

U^,^\r)^Cp'^, (14) 

Wr(r)-OAro(r). (15) 

The terms A, C, and D are scalars. Hence, their form is natural for all shapes of 
the nucleus. The term B is valid only in case of spherical symmetry. It is inspired 
by the Ma and Wambach parameterization of lEM. 

The terms A and B depend on the first derivative of the particle density. This 
derivative is peaked at the surface, so these terms are expected to yield the desired 
profile of the lEM. Since the gradient square is positive and the radial derivative 
itself is negative, an upward-pointing peak will require negative and positive values 
of the coupling constants in cases A and B, respectively. 

The term D depends on the second derivative. Hence, it should produce two 
peaks of opposite signs at the borders of the surface, where the curvature of the 
density profile is largest. 



3. Radial profiles of the isoscalar effective mass 

Correct description of the s. p. -level density near the Fermi surface requires the 
mean value of the lEM to be close to unity. In order to fulfill this condition we 
impose the following constraint on the radial lEM profile, 

(16) 

J A m ^ ' 

where A is the mass number of the considered nucleus. 

We examine each of the four new terms, A-D, separately, that is, by switching 
the remaining three off. In each case, we vary the concerned new coupling constant 
by hand, and for each value thereof, we readjust the coupling constant to fulfill 
the condition (fT6|l . This means that the excess of the lEM produced by the peak is 
compensated by decreasing the lEM inside the nucleus. 

The remaining coupling constants of the functional ([2]) are kept intact at values 
of the SkXc Skyrme parameterization^l^l This specific parameterization was chosen 
because of the condition (fTHll . Indeed, in order to fulfill Eq. (fT6|) it is reasonable 
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Fig. 1: Radial dependence of the JEM, mQ{r)/m, in ^°Ca for the four variants of 
the calculations. The solid and dashed lines show the lEM profiles for the limiting 
achievable values of the corresponding coupling constants. See text for details. 

to choose a parameterization with the lEM close to unity, which is a case for the 
SkXc. In addition, this force was fitted with particular attention payed to the s.p. 
spectra. 

It turned out that the possible values of the new coupling constants are limited. 
Extending the new coupling constants beyond certain limits leads to strong oscilla- 
tions of the density what causes that the total energy diverges. The radial profiles 
of the lEM obtained in ""^"Ca for the limiting (both negative and positive) values 
of the new coupling constants are shown in Fig. [1] The limits corresponding to 
downward-pointing peaks in the lEM are, most likely, of limited physical interest. 

It is clearly seen that in the variants A and B of our model the desired surface- 
peaked lEM profile is obtained. In the variant C a well pronounced peak appears 
as well, but it is shifted toward the center of the nucleus. The variant D yields only 
small fluctuations of the lEM value, and it does not seem to be of much relevance. 

4. Influence of the new terms on single-particle energies 

In our study, we are mostly interested in the influence of the new terms on the 
s.p. spectra, in particular on the spin-orbit (SO) coupling. For each pair of the 
SO partners, with the same n, I, and j = I ± ^, we define the SO splitting as the 
difference of their energies, and the SO centroid as the arithmetic average thereof. 
Figure [2] shows, separately for each variant of our model, the SO splittings and 
centroids as a function of new coupling constant for the entire accessible range of 
each coupling constant. 

The new terms influence both the SO splittings and centroids quite strongly in 
all variants of our model except the variant D. In the latter case, the modifications 
of the lEM's radial profile are apparently too weak to produce any clear trend. In 
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Fig. 2: Spin-orbit splittings (upper panels) and centroids (lower panels) as a function 
of C^37 (a) and (b), C^' (c) and C^p (d) in ^OQa (left), ^^Ni (middle) and 

i°°Sn (right) nuclei. 



our calculations, the term C affects the s.p. levels most strongly. However, it also 
shifts the binding energies of the considered nuclei by more than 50% when going 

2 

from Cq = to the limiting values. It is, therefore, clear that in order to remain 
within a physically acceptable area, it would be necessary to refit all the coupling 
constants of the functional at least to masses and radii of the three considered 
nuclei. This is, however, beyond the scope of the present work. Such a problem is 
of slightly lesser importance for the terms A and B, since their inclusion does not 
change the masses by more than 10%. 

As expected, the terms A and B give a surface-peaked lEM profile for negative 
and positive values of their coupling constants, respectively. Taking this into ac- 
count, one can see from Figs. [2] (a) and (b) that the impact of the peak on the SO 
splittings is opposite in these two cases, although the plots look similarly. For the 
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term A, the emergence of the peak leads to an enhancement of the SO sphtting of 
the If level in *'^Ca. This trend is undesired since all the existing Skyrme parar ne- 
terizations overestimate this splitting already without additional terms, see Ref. 
and articles quoted therein. On the other hand, the onset of the peak caused by 
the term B quenches the If splitting in ^''Ca, driving it closer to the experimental 
value. 

5. Summary and conclusions 

We presented a way to extend the standard Skyrme EDF by introducing new terms: 
t(Vp)^, r^, r^, and rAp. The rationale behind is to modify the radial profile of 
the lEM toward a surface-peaked geometry. In this way one can effectively include, 
within the nuclear EDF theory, coupling to surface vibrations and restitute correct 
level density at the Fermi surfacel^^. Our study shows, in particular, that the term 
gives the desired surface-peaked lEM and that it seems to be most promising 
as far as reproducing the experimental s.p. spectra is concerned. However, this term 
is not a scalar, and must be therefore generalized in order to be applicable to non- 
spherical shapes. A more detailed analysis of the variants A and B of our model is 
underway. 
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